function [ relatedGenes ] = getGenesForTrait( startSite, endSite, geneNames, geneChrs, geneStarts, geneEnds )
%UNTITLED2 Summary of this function goes here
%   Detailed explanation goes here
    
    numGenes = size(geneChrs, 1);
    
    relatedGeneIdx = 0;
    for i = 1:numGenes
        if (geneStarts(i, 1) > startSite || geneEnds(i, 1) < endSite)
            relatedGeneIdx = relatedGeneIdx + 1;
            relatedGenes{relatedGeneIdx, 1} = geneNames{i, 1};
        end
    end

end

